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Abstract 



The paper deals with reaction-diffusion equations involving a hysteretic discontinuity in the source 
term, which is defined at each spatial point. Such problems describe biological processes and chemical 
| reactions in which diffusive and nondiffusive substances interact according to hysteresis law. Under 

the assumption that the initial data are spatially transverse, we prove a theorem on the uniqueness 
of solutions. The theorem covers the case of non-Lipschitz hysteresis branches arising in the theory of 
slow-fast systems. 
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1 Introduction 

We consider reaction-diffusion equations with right-hand sides involving a discontinuous hysteresis defined at 
each spatial point. Such problems describe biological processes and chemical reactions in which diffusive and 
nondiffusive substances interact according to hysteresis law. As a result, various spatial and spatio-temporal 
patterns may appear (see, e.g., 015]). 

First rigorous results about the existence of solutions of parabolic equations with hysteresis in the source 
term have been obtained in [UEB 16 f° r multi-valued hysteresis. Formal asymptotic expansions of solutions 
were recently obtained for some special case in [§]. However, the uniqueness of solutions and their continuous 
dependence on initial data as well as a thorough analysis of pattern formation remained open questions. 

In [3J, a new approach was suggested. It allowed us to find a broad class of initial data (transverse 
functions, see Sec. H]) for which a solution exists and, if unique, continuously depends on initial data. The 
approach is based on tracking the so-called free boundary which defines the hysteresis topology. The hys- 
teresis topology is, in the simplest case, related to the structure of subdomains in space where hysteresis 
takes the same value. The main advantage of this approach is that, tracking the hysteretic free boundary, 
one gets necessary information on the precise form of emerging spatio-temporal patterns. 

The transversality assumption roughly speaking means that the initial function has a nonvanishing deriva- 
tive on the boundary between the above-mentioned subdomains at the initial moment (see Condition 12.21) . 
Under this assumption, we prove in the present paper that the solution for the reaction-diffusion equation 
with discontinuous spatially distributed hysteresis is unique for some positive time. Combining this result 
with the results in [5J, we obtain that the solution can be uniquely continued as long as it remains transverse. 

We consider a one-dimensional domain and, for the clarity of exposition, we concentrate on a scalar 
reaction-diffusion equation. The paper is organized as follows. In Sec. [3J we define functional spaces, 
introduce spatially distributed hysteresis and set the prototype problem. In the end of Sec. [5J we formulate 
the main result of the paper: Theorem l2.Il on the uniqueness of transverse solutions. 
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Section [3] is devoted to the proof of Theorem 12.11 

Interestingly, the uniqueness of solutions holds for some classes non-Lipschitz hysteresis branches, too. 
In particular, we prove the uniqueness for hysteresis branches arising in bistable slow-fast reaction-diffusion 
systems (see e.g., |13j and the references therein). In Appendix [K\ we briefly describe such systems and show 
that the arising hysteresis branches satisfy assumptions of our main theorem. 



2 Setting of the problem 

2.1 Functional spaces 

We denote by L q = L q (0, 1), q > 1, the standard Lebesgue space and by W q = W q (0, 1) with natural I the 
standard Sobolev space. For a noninteger I > 0, denote by W q = W q (0, 1) the Sobolev space with the norm 



, r 1 , r 1 i«(ra>(z)-i;(ro%)|9 



1/9 

\ x - y|l+9('-W) 



where [I] is the integer part of I. 

Let Qt = (0, 1) x (0, T). We introduce the Holder space C' f (Q T ), < 7 < 1, and the anisotropic Sobolev 
space W q ' (Qt) with the norm 



\\ u Ww^(Q T ) = / \H;t)\\w*dt+ I ||ttt(-,t)||i„tft 



1/'/ 



Throughout the paper, we fix q and 7 such that 

q>3, 0< 7 <l-3/<7. 

Then u 7 u x £ C 7 (Q T ) whenever u £ W q 1 (Qr) (see Lemma 3.3 in [12j Chap. 2]). 

In what follows, we will consider solutions of parabolic problems in the space W q ,1 (QT)- To define the 
space of initial data, we will use the fact that if u £ W^' (Qt), then the trace u|t=t is well defined and 
belongs to W q 2/>q for all to £ [0, T] (see Lemma 2.4 in [12l Chap. 2]). Moreover, one can define the space 

^ q as the subspace of functions from W q 2 ^ q with the zero Neumann boundary conditions. 

2.2 Hysteresis 

In this section, we introduce a hysteresis operator defined for functions of time variable t. Then we extend 
the definition to a spatially distributed hysteresis acting on a space of functions of time variable t and space 
variable x. 

We fix two numbers a and j3 such that a < /3. The numbers a and /3 will play a role of thresholds for 
the hysteresis operator. Next, we introduce continuous functions (hysteresis branches) 

Hi : (-00, P] i-> K, H 2 : [a, 00) ^ E. 

We assume throughout that the following condition holds. 

Condition 2.1. There is a number a £ [0, 1) such that, for any U > 0, there exists M = M(U) > with 
the properties 

M 

\Hi(u) - ffx(u)l < _ \u - u\, Vu,u£[-U,p), (2.1) 

M 

\H 2 (u)-H 2 (u)\<- — — Mu,u£(a,U]- (2-2) 

(u — a) a + (u — a) a 
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Remark 2.1. 1. Any locally Lipschitz continuous functions Hi(u) and H2{u) satisfy Condition 12.11 
Moreover, this condition covers the important case where the hysteresis branches H±(u) and H2(u) are 
the stable parts of the curve g(u, v) = in the slow-fast system (|A.1[) . They are not Lipschitz contin- 
uous near the points u = /3 and a, respectively. But they still satisfy Condition 12.11 (see Appendix [X] 
for details). 

2. On the other hand, any Hi(u) and H2(u) satisfying Condition 12. II are locally Holder continuous with 
exponent 1 — a on (-co, 0\ and [a, oo), respectively. Furthermore, H±(u) and H2(u) are locally Lipschitz 
on the open intervals (— oo, /3) and (a, oo), respectively. Therefore, they satisfy the assumptions in [3J, 
where existence of solutions and their continuous dependence on initial data were proved. 

We fix T > and denote by C r [0,T) the linear space of functions which are continuous on the right 
in [0,T). For any £o € {1,2} (initial configuration of hysteresis) and g G C[0,T] (input), we introduce the 
configuration function 

C : {1,2} x C[0,T] C r [0,T), CW = C(Co,ff)(t) 
as follows. Let X t = {f G (0,i\ : g(t') = a or /?}. Then 

!1 if ,g(0) < a, 
2 if .g(0) > /3, 
Co if 5(0) G (a,p) 

and for t G (0, T] 

ifX t = 0, 

if X t ^ and g(maxA t ) = a, 
if Xt 7^ and g(maxA t ) = 

Now we introduce the hysteresis operator (cf. [TTin5] ) 

H : {1,2} x C[0,T] -> C r [0,T) 

by the following rule. For any initial configuration Co G {1,2} and input g G C[0,T], the function H(Co,g) '■ 
[0,T] -> R (output) is given by 

H(( ,g)(t)=H m (g(t)), 

where ((t) is the configuration function defined above (see Fig. 12. 1| . Note that this hysteresis operator is 
typically not continuous. 



at) 





Figure 2.1: The hysteresis operator % 

Now we introduce a spatially distributed hysteresis. Assume that the initial configuration and the input 
function depend on spatial variable x G [0, 1]. Denote them by £,o( x ) an d w(a;, t), where 



£ : [0,1] H>{1,2}, u:[0,1]x[0,T]h> 
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Let u(x, •) € C[0, T], Denote = 0). We say that a function <p(x) and a hysteresis configuration 
£o{x) are consistent if, for any x £ [0, 1], 

[{1} X<p{x)>0, 
({1,2} iiip(x)e(a,ft). 

Then we can define the function 

v(x,t)=n{£ (x),u(x,-))(t), (2.3) 
which is called spatially distributed hysteresis. 

2.3 Reaction-diffusion equations with hysteresis 

The main object of this paper is the initial boundary-value problem for the reaction-diffusion equation 

u t = u xx + v, (x,t)€Qr, (2.4) 
u x \ x =o = u x \x=i = 0, (2.5) 
u\t=Q = (p(x), 16(0,1), (2.6) 

where v = v(x, t) represents the spatially distributed hysteresis given by (|2.3[) . 
Remark 2.2. The more general equation 

u t = u xx + f(u,v) 

with locally Lipschitz continuous right-hand side / can be reduced to Eq. (|2.4[) . Indeed, it suffices to replace 
the hysteresis branches Hj(u) in the definition of hysteresis H by Fj(u) = f(u,Hj(u)), j — 1,2. One can 
check that the functions Fj (u) also satisfy Condition 12.11 

The general assumption on the initial data <p(x) under which the uniqueness result holds is that <p{x) is 
transverse with respect to the initial configuration £o(x). 

Definition 2.1. We say that a function <p e C 1 ^,!] is transverse (with respect to a spatial configuration 
£,o(x)) if it is consistent with £o(x) and the following holds: 

1. if ip(x) = a and ip'(x) — for some x € [0,1], then £o(x) = 1 in a neighborhood of x; 

2. if ip(x) — (3 and <p'(x) = for some 16 [0,1], then £,q(x) — 2 in a neighborhood of x. 

For the clarity of exposition, we now restrict ourselves to the following prototype situation (see Fig. 12.21) . 
Fix some a £ (0, 1). 

Condition 2.2. 1. Forae (0,1), one has 

to(x) = < 

12, x > a. 

2. The equation ip(x) = ft on the interval [0,a] has no roots. 

3. The equation (p(x) — a on the interval [a, 1] has the unique root x = a. 

4. <p'(a) > 0. 
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Figure 2.3: Initial data satisfying the assump- 
Figure 2.2: Initial data satisfying Condition O tions of Remark [231 



It follows from this condition that the hysteresis (|2.3j) at the initial moment is given by 

, q = iH 1 (ip(x)), x < a, 
\H 2 (tp(x)), x>a. 

We give a definition of a solution of problem (|2.4p ~ (|2.6[) . assuming that ip £ ^ q . 

Definition 2.2. A function u(x,t) is called a solution of problem (|2.4I) - (|2.6[) (in Qt) if u £ W%' (Qt)> 
v(x,t) is measurable, u and v satisfy equation (|2.4I) for a.e. (x, i) € Qt, and conditions (|2.5[) and (|2.6I) arc 
satisfied in the sense of traces. 

It follows from this definition that any solution u(x, t) belongs to C(Q T ). Therefore, the function v(x,t) 
is well defined by (12.31) and belongs to L oq (Qt)- 

Remark 2.3. An initial function may also satisfy (p(a) > a. Such a function is also consistent with £q(x) 
from Condition l2.2[ provided that the equation tp(x) = (3 has no roots on the interval [0, a] and the equation 
tp(x) = a has no roots on the interval [a, 1] (see Fig. I2.3[) . However, this situation is much simpler. Indeed, 
one can show (see [3]) that, for any solution of problem (|2.4[) - (|2.6|) with such an initial function, the hysteresis 
1-L(^q(x),u(x, -))(t) does not depend on t on a sufficiently small time interval. Therefore, in this case, we 
leave details to the reader and rather concentrate on Condition 12.21 

2.4 Main result 

In what follows, we always assume that Conditions 12.11 and 12.21 hold. 

In [3], it was proved that any solution of problem ()2.4p - (|2.6j) remains transverse on some time interval 
[0, T] in the sense of the following definition. 

Definition 2.3. A function u £ C 1,0 (Q T ) is transverse on [0,T] (with respect to a spatial configuration 
£(x, t)) if, for every fixed t £ [0, T], the function u(-,t) is transverse with respect to the spatial configuration 

£(•,*)■ 

Furthermore, problem (|2.4|) - (|2.6|) has at least one transverse solution on the above interval [0, T] and any 
such solution can be extended to its maximal interval of existence on which it remains transverse. 

In this paper, we show that there exists no more than one transverse solution of problem (|2.4[) - (|2.6|) . We 
formulate the main result as follows. 

Theorem 2.1. Let Conditions 12.11 and 12.21 hold, and let q > 3. Assume that u,u £ Wg ' (Qt ) are two 
transverse solutions of problem (|2.4[) - (12.61) in Qt for some Tq. Then u = u. 
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3 Proof of Theorem O 



Taking into account Rcmark l2.31 it suffices to prove the uniqueness on a sufficiently small time interval [0, T], 
T < T . 

We denote _ 

_ fx (a) 

By Condition^ f > 0. 

Lemma 3.1. Let u G Wg ' (Qt ) be a solution of problem (j2.4p - (|2.6[) for some Tq > 0. Then there exist 
T € (0,To) and 5 > suc/i f/iaf the following hold on the interval t G [0,T]: 

1. u x (x, t) > f for all x € [a — <5, a + 5]: 

2. the equation u(x, t) = a on the interval x G [a — S, 1] has a unique root x = a(t); 

3. a(t) is continuous; 
4- the function 

b(t) = max a(s) (3.1) 
«e[o,t] 

satisfies b(t) G [a, a + 8]; 

5. the equation u(x,t) = ft on the interval x G [0, b(t)] has no roots. 

Proof. The assertions of the lemma follow from the fact that it, u x G C~'{Q T ) for all sufficiently small T (with 
the norms in C 7 (Q T ) bounded uniformly with respect to small T), u(x,0) = f{x), and f(x) is transverse 
(see Condition |2~2")) . □ 

By possibly decreasing T and 5, we see that Lemma 13.11 holds for u with some functions a(t) and 

b(t) = max a(s) instead of a(t) and bit). 
«e[o,t] 

Now the key observation is that the hysteresis H. acting on transverse functions it and u on the time 
interval [0, T] can be represented in terms of the free boundaries b(t) and b(t) as follows: 



H(£o(x),u{x,-)(t) 
H(£ (x),u(x,-)(t) 



Hi(u(x,t)), 0<x<b(t), 

H 2 (u(x,t)), b(t)<x<l, 

: . (3-2) 

Hi(u(x,t)), 0<x<b{t), 

H 2 (u(x,t)), b(t)<x<l. 



The following lemma allows us to estimate the distance between the free boundaries b(t) and b(t) in terms 
of the difference between u and u. 

Lemma 3.2. Let T be the number from Lemma |3. 11 Then 



\b- b\\ C [o,T\ < =11" - u 



1 m 

11 " C(Q T )- 



f 

Proof. It follows from the definition of b(t) and b(t) that 

\\b - b\\c[o,T] < \\a - a||c[o,T]- 
On the other hand, using Lemma 13. 11 we obtain for any t G [0,T] 

\a(t)-a(t)\ < ku(a(t),t) - u(a(t),t)\ < ^\\u - u\\ c( q t) , 
which completes the proof. □ 



G 



Now we can prove Theorem 12.11 

1. Denote w = u — u. The function w satisfies the linear parabolic equation 

m = w xx + h(x, t), (x, t) £ Q T , (3.3) 

where h(x,t) = 'H(^q(x),u(x, •)) — H(t;o(x),u(x, •)), and the zero boundary and initial conditions. Obviously, 
h £ Loo{Qt), and the function w can be represented via the Green function G{x, y, t, s) of the heat equation 
with the Neumann boundary conditions: 

t l 

w(x,t) = [ [ G(x,y,t, s)h(y, s) dyds. 

Therefore, using the estimate 

\G(x,y,t,s)\ < -J^ e -(-*) 2 /(4(*-*)) < < s < t, 

with fci > not depending on (x,t) g Qt and T > (see, e.g., 0), we obtain 

t l 

\w(x,t)\ < fci / -=L= / h(y,s)dy. (3.4) 



o o 





2. Let us estimate the interior integral in (|3.4I) for a fixed s. We assume that b(s) < b(s) (the case 
b(s) > b(s) is treated analogously). Then, due to (|3.2[) . 

fffiCtO-fliCu), 0<y<6( S ), 
%,s) = < ff 2 («) - 6(s) < y < b(a), 

[h 2 {u)-H 2 {u), b(s)<y<l. 

2.1. Assertion 5 in Lemma 13-11 implies that 

u(y,s) < f3, u{y,s) < (3 

on the closed set {(y,s) : y £ [0,b(s)}, s £ [0,T]}. Hence, the values f3 — u(y,s) and f3 — u(y,s) are separated 
from 0. Therefore, using Condition 12.11 we obtain 

b(s) b(s) 

M 



J W - u{y , s)r + {f 3-u( y , s ))>^ s) -^ s)ldy 



b(s) 

< k 2 J \u(y,s)-u(y,s)\dy <k 2 \\w\\ c( Q T) , (3.5) 
o 

where k 2 > and the constants ^3,^4,^5 > below do not depend on s £ [0, T]. 
2.2. By the boundedness of Hi(u) and H 2 (u) for (y,s) £ Q T , we have 

b(s) b(s) 

\h(y> s )\dy < h J dy < ka\\b - b\\c[o,T\- 

b(s) b(s) 
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Applying Lemma I3T21 yields 

$0) 

£3 



&0) 



\h(y,s)\dy<^\\w\\ c( Q T) . (3.6) 



2.3. Let 5 be the number from Lemma [3.1[ and let b(s) < y < a + 5. Then, using assertion 1 in Lemma |3. II 
and the mean- value theorem, we have (see Fig. I3.1[) 



u(y,s 




Figure 3.1: The mean-value theorem for u(y,s): tancj > ip 



\u(y, s)-a\= u{y, s) - u(a(s), s) > (y - a(s))(p > (y - b(s))ip. 

Similarly, 

\u(y,s)-a\ > {y~b{s))Tp. 
Taking into account these two inequalities and using Condition 12.11 we obtain 

l+S \ h (y, S )\ d y < M- a+ f 5 Hy,s)-u(y,s)\ 

' ^ n y ~ 2 J {y-b{s)Y 

b(s) b(s) 



rii 11 a+5 
M W W \\c(Q T ) f I 



(3.7) 



2 

H 

2.4. Finally, assertions 1 and 2 in Lemma |3~T1 imply that 



(«-&(«)) 



-dy< fc 4 ||Hlc(Q T )- 



u(y,s) > a, u(y,s) > a 

on the closed set [a + 6,1} x [0,T]. Hence, the values u(y,s) — a and u(y,s) — a are separated from 0. 
Therefore, due to Condition I2JJ 



(3.8) 



1 1 

f M 

\h{y,s)\dy< / -— ; 7 . jy, r —\u(y,s)-u{y,s)\dy 

_J {u{y, s) - ay + {u{y, s) - ay 

a-\-8 a+5 

1 

< h J \u(y,s)-u(y,s)\dy <k 5 \\w\\ c( Q T y 

a+5 

3. Combining estimate (|3.4p with inequalities (|3.5p - (|3.8p . we obtain 



1 ds 



\w(x,t)\ < k 6 \\w\\ c( Q T) J -^==2k 6 T 1 / 2 \\w\\ c{ Q T} 
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Since k$ does not depend on x and t, the latter inequality is equivalent to ||Hlc(Q T ) — ^^ 1 ^ 2 \\ w \\c(Q T )' 
Since ke does not depend on (small) T either, it follows that w = provided that T > is small enough. □ 

A Connection with slow-fast systems 

It is known [3UTU] that hysteresis may approximate solutions of ordinary differential equations with a small 
parameter e > and a "bistable" right-hand side. Combining such an ordinary differential equation with 
the reaction-diffusion equation (|2.4[) . one obtains a slow- fast system of the form 



u t = u xx + f(u,v), 
ev t = g(u,v), 



(A.I) 



where e > 0, f(u, v) and g(u, v) are smooth functions of class C°°, and g(u, v) satisfies the following condition 
(see Fig. IA.1I) . 



9 fa, v) 








Figure A.l: The nullcline of g(u,v) 



Condition A.l. 1. There are two numbers a < ft such that the equation g(u,v) = has a unique root 
v = H\(u) for u < a and v = Hq.(u) for u > /3; three distinct roots v = H\{u), v — H2(u), and 
v = Hq(u) for u £ (a, ft); two distinct roots v = H\(a) and v — H2(a) = Ho(a); two distinct roots 
v = Hxifi) = H (P) and v = H 2 (/3). 

2. the functions H\{u) and H2(u) are locally Lipschitz continuous for u < ft and u > a, respectively. 

3. Let A = (a,H 2 (a)) and B = (ft, HA ft)). Then -i^S. + o and + 0. 

ou ou 



8g(A) _ _ d^gjA) _ Q d n g(A ) 
dv dv 71 " 1 dv r 



2. — = • • • = — z — = 0, — ^ for some even n > 2 and the same holds at the point B. 



5. g(u,v) > (< 0) if a point (u,v) lies to the left (right) from the curve g(u,v) — on the plane (u,v). 

The simplest example of such a bistable nonlinearity typical, e.g., for the FitzHugh-Nagumo systems is 
given by 

g(u, v) — u + v — v 3 . 

As we have mentioned before, the natural conjecture is that the solutions of the slow-fast system (|A.1|) 
approximate as e — > the solution of Eq. (|2.4|) with hysteresis H defined by the curves H±(u) and H2(u) 
from Condition IA.ll We refer to [13] , where singular limit analysis of traveling waves in bistable systems is 
done, and to [2lll4j. where equations of the form u t — A$(it) with a cubic nonlinearity $ are studied. 
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The proof of the above conjecture in the general situation is an open question, which is beyond the scope 
of the present paper. However, it is clear that the proof would deal with hysteresis defined by non-Lipschitz 
curves. The following result assures that the non-Lipschitz curves Hi(u) and H2(u) from Condition IA.1I 
satisfy Condition 12.11 and thus fit into our theory. 



Lemma A.l. Let g(u,v) satisfy Condition \A.1\ Then the functions H\{u) and H2{u) from Condition \A.\\ 
satisfy Condition \'2.1\ with a — (n — l)/n. 

Proof. 1. Without loss of generality, we assume that A = and prove the lemma for the function H^iu). 
Due to item 2 in Condition IA.1[ it suffices to prove inequality (|2.2|) for u in a small neighborhood of a = 0. 

Since — - — 7^ 0. it follows from the implicit function theorem that there is a unique function G2M 
du 

defined in a small neighborhood of such that g{Gi{v), v) = and 

G 2 (0) = G' 2 (0) = ■■• = G 2 ™~ 1) (0) = 0, G 2 n) (0)>0. (A.2) 

Obviously, G^v) and H^iu) are inverse to each other for all negative v and positive u close to 0. It is 
convenient to introduce the function G(w) — G 2 (— 1«) defined for small positive w. Due to (|A.2|) . it satisfies 

G(0) = G'(0) = ...G ( ™- 1) (0) = 0, G (n) (0)>0 (A.3) 

since n is even. Now if we denote w = —H 2 (u) (> 0), then u = G2(H2(u)) = G 2 (— u>) = G(w). Therefore, 
the inequality for H2 in Condition (|2.1[) . which we have to prove, is equivalent to the following: 



m - m > M 



([GW]^ + [GH]^) (A.4) 



for all < w < w < eq, where Eq > is sufficiently small and M > does not depend on w, w. 
2. Expanding G(w) by the Taylor formula about w = w, we have 

m^m = g'{w) + g ^m (w _ ^ + _ ^ (A . 5) 

w — w z — ' fc! TV. 

k=2 

where < £ = £(w, w) < Eq. 

Further, we expand G^ k '(w), k = 1, . . . , n — 1 by the Taylor formula about the origin, using (|A.3[) : 

[n — k)\ 

where < — £k(u>) < £o- Substituting these expressions into (IA.5I) yields 

w — w (n — 1)1 n\ , , , 



2(n — lj! \ n 

where we assume e so small that G^(6c) > ^ and G< n >(£) > 

It is easy to show that 

w™- 1 + - w)"- 1 > Mi (to"- 1 + u-™- 1 ) 
n 

for some Mi > not depending on w and w, < w < w. Hence, inequality ()A.6|) implies 

° H - > M 2 + ^) , (A.7) 

w — w 
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where M2 > does not depend on w and w, < w < w < e. Combining (IA.7I) and the relations 

GH = < 

n! n! 

with < C = CM < £0 yields pO|) . □ 
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